Metal nanoparticles with sharp corners: Universal properties of plasmon resonances 
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We predict the simultaneous occurrence of two fundamental phenomena for metal nanoparticles 
possessing sharp corners: First, the main plasmonic dipolar mode experiences strong red shift with 
decreasing corner curvature radius; its resonant frequency is controlled by the apex angle of the 
corner and the normalized (to the particle size) corner curvature. Second, the split-off plasmonic 
mode experiences strong localization at the corners. Altogether, this paves the way for tailoring of 
metal nano-structures providing wavelength-selective excitation of localized plasmons and a strong 
near-field enhancement of linear and nonlinear optical phenomena. 
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Plasmon excitations of metal nanoparticles, including 
2D particles (nanowires) , is a vast and hot research area. 
Potential applications of nano-plasmons span from nano- 
lasers (spasers) [l|-[3| to bio-sensors and sensors of sin- 
gle atoms and molecules 0, HJ. It is highly important 
to apprehend the possibilities of tailoring of plasmonic 
resonances and the extends of concentration of the light 
energy deeply on the sub-wavelength scale. 

Analytical plasmonic solutions to Maxwell's equations 
are available only for metal nanoparticles of simplest 2D 
and 3D shapes - like circular cylinder, sphere, and ellip- 
soid [1, 0] • They include the optical permittivity of the 
metal £m(w) = e^fw) + ie'^co), such that e' M < and 
ejjj <C \s'm\i an d give often strongly degenerated plas- 
monic resonances. In particular, all plasmonic excita- 
tions of a circular cylinder correspond to a single eigen- 
frequency given by ey(w) = — 1. 

Using direct numerical methods, serious efforts have 
been undertaken to ascertain the phenomenology of the 
plasmonic response of 2D and 3D particles of more com- 
plicated shapes @-ll|. The general outcome is that the 



deviations from the most symmetric shapes result in en- 
riching plasmon spectrum and in the appearance of new 
red-shifted resonances. In the 2D case, the strongest 
shape effects were found for triangular cross-sections. 

Regardless of the plasmonic effects, the presence of 
sharp tips and corners (metal or dielectric) is known 
to lead to the corner singularities, i.e., to a strong lo- 
cal enhancement of electromagnetic fields [ij], E3]- At 
the same time, the impact of tips and corners on the 
plasmonic properties of nanoparticles remains in essence 
unexplored. In numerical simulations, the sharp surface 
features are typically rounded to avoid numerical insta- 
bilities. Non-rounded 7r/2-corners lead to divergence of 
numerical methods for £m — > —3 II 41 . Mathematical 
foundations of the plasmonic theory |l5l . Il6j refer to suf- 
ficiently smooth (Lyapunov) surfaces. A recent heuristic 
attempt to deal with plasmonic properties of nanopar- 
ticles possessing perfectly sharp features 117] is strongly 
unconvincing. 

Experimentally, the presence of sharp features of metal 



nanoparticles is far from exotic. Modern chemical meth- 
ods for production such particles are able to provide al- 
most atomically sharp edges in the crystallographic direc- 
tions [l8l - |20j . Synthesis of metal nanostars has allowed 
recently to stron gly enhance the near fields at the plas- 
monic resonance [21| . 

In this Letter, we predict and analyze a strong im- 
pact of corners of metal nanoparticles on fundamental 
plasmonic properties. Both the spectrum and the spatial 
structure of the plasmonic modes are strongly affected by 
the corners showing simple and universal features. 

We employ the most efficient approach to the descrip- 
tion of plasmonic resonances - the method of integral 
equations for the amplitude of the surface-charge density 
cr(r) [22, 23 1 . Applicable in the quasi-static approxima- 
tion, it treats the plasmonic eigenproblem as a geometric 
issue: The real resonant values of the optical permittiv- 
ity Sj depend only of the shape of the particle, while the 
corresponding eigenfunctions <7j(r) are scaling invariant. 
The knowledge of the eigenmodes allows to determine the 
width of the resonance, the polarizability of the particle, 
near fields, etc., using a perturbation routine. Reduction 
of the dimension of the problem to be solved strongly 
enhances the capability of numerical methods. 

In the 2D case, the eigenproblem to be solved reads: 
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where the line L is the boundary of the particle, dl is the 
length element along L and Kj — [sj + 1) / (ej — 1) is the 
eigenvalue. The kernel K is real and given by 



K(t,t') 



n • (r 



7T (r 



(2) 



where n = n(r) is the unit vector of the external normal. 
Generally, Eq. (JXJ) gives an infinite sequence of the modes. 

The integral eigenproblem for Ej and cr,-(r) possesses 
important general features [23}. The total modal charge 
is zero, J L aj(r)dl — 0. Since the kernel is not sym- 
metric, K(r,r') ^ K(r',r), the eigenproblem is not 
Hermitian. Thus, one has to employ additionally the 
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eigenfunctions Tj(r) of the adjoint problem with the 
transposed kernel K (r', r) and the orthogonality relation 
f L <jj(r) Tj'(r) dl = for j ^ j' . In the 2D case, the 
eigenvalues appear in twin pairs e j , £ J 1 . 

We specify the boundary line L by the polar-angle de- 
pendence of the radius, \r\(<p) = r(ip). To investigate the 
impact of corners, we employed different parametriza- 
tions of r(ip) with variable corner-curvature radius p c and 
apex angle 9 a < n/2. The simplest one is given by 



solid line for 9 a — tt/2 tend to 1: This is the limit of the 



circular cross-section with e, 



T. We have made sure 
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where r = r(0) is the large half-diameter, a = ctg(# a /2), 
and 1 < p < oo. It gives a smoothed rhombus with the 
sharp apex angle 9 a and the normalized corner curvature 
k c = r /p c = ctg 2 (6> a /2)(p - 1 growing linearly 

with p for p> 1, see also Fig. [TJl The case p — 1 cor- 
responds to an ellipse with the axes ratio of a. Another 
parametrization, illustrated by Fig. [T|d, is a smoothed 
equilateral triangle (8 a — 7r/3) transferring to the cir- 
cle for p c — tq. The kernel K(r, r') can be expressed 
by r(tp), r{(p') 1 and dr/dip; it peaks at the corners with 
K(r c ,r c ) = l/2np c . 




FIG. 1. The boundary line L for a smoothed rhombus trans- 
ferring to an ellipse (a) and a smoothed equilateral triangle 
transferring to a circle (b). 

Equation (TTJ) with the described boundary lines L was 
solved numerically using crowding of the points near the 
corners. It was made sure that further increasing of the 
number of points does not influence the results. All gen- 
eral properties of £j and aj (r) were fulfilled. 

The solid lines in Fig. [5] show the dependences £i(k c ) 
for the lowest dipolar branch and several values of the 
apex angle 9 a of the rhombus on a semi-logarithmic scale. 
The starting values of E\ correspond to the known plas- 
monic solutions for the elliptic cross-section. Increasing 
corner curvature k c results in strong decrease of the res- 
onant permittivities [i.e., in strong red shifts of the res- 
onant frequencies in accordance with the Drude-like de- 
pendences £m(w)] for all values of 9 a - This decrease per- 
sists for k c > 10, when the shape is already settled down 
and increasing sharpness of the corners remains the only 
variable feature. The dotted line is plotted for the trian- 
gular parametrization, a — 7r/3. When k c increases, it 
approaches quickly the solid line corresponding to rhom- 
bus with the same 9 a . For k c — > 1, the dotted line and the 



also that different types of smoothing give essentially the 
same results for k c 3> I. 

The dipolar branches £j(k c ) with j > 2 are well 
separated from the main branch £i(k c ) for k c ^> 1, 
£2,3 1 ~ + This makes possible selective excita- 
tion of the main branch. At the same time, the vicinity 
of the point e = — 1 is always filled up with the spectrum. 
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Normalized corner curvature, s: = r I p c 

FIG. 2. The lowest dipolar branch ei(k c ) for the smoothed 
rhombus with the apex angle 9 a = 7r/2, 7t/3, 7r/4, tt/5, and 
7r/6 (the solid lines) and a smoothed equilateral triangle with 
8 a — 7r/3 (the dotted line). The horizontal dashed lines show 
the corresponding limiting values 1 — 2-k jQ a . 

The horizontal dashed lines in Fig. [^indicate the crit- 
ical values e c = 1 — 2ir/9 a which correspond to non- 
integrable field singularities for a perfect single corner 
with the apex angle 9 a [Hj]- Obviously, the branch Si{n c ) 
cannot cross the corresponding horizontal line. At the 
same time, the tendency of approaching the critical val- 
ues is clearly seen, especially for not very small values of 
9 a . This approaching strongly slows down with increas- 
ing k c so that the distances Si(K c ,9 a ) — e c {9 a ) remain 
larger than (0.2 — 0.8) even for n c w 10 5 , i.e., for non- 
realistic sub-atomically sharp corners. This shows that 
the perfect-corner limit is of no practical importance. 

Next, we analyze the impact of varying corner curva- 
ture on the spatial structure of the eigenfunction <7i(r), 
representing the surface-charge density for the lowest 
dipolar mode. When using the polar angle ip, it is 
practical to employ the angular charge density (Jj(p) — 
tjj[r(tp)]dl/d(p, such that <7j(ip)dip gives the charge differ- 
ential. Owing to non-hermitian nature of the plasmonic 
eigenproblem, normalization of <Jj(tp) is a matter of con- 
venience. Our first choice is CTj(O) = 1. 

Consider the representative case of rhombus with 9 a = 
7r/3; it corresponds to non-degenerate branch Sj(K c ). 
The solid lines in Fig. [3^ show the dependence <Ji(p) 
for three values of k c . All necessary properties of spatial 
symmetry are fulfilled. The sharp corners (p = and tt 
are charged, while the obtuse corners ip = ±tt/2 are un- 



3 




cp In 



FIG. 3. (a) The eigenfunction <Ti(<p) for the rhombus with 
9 a ~ 7r/3 for n c = 2, 4, and 16. (b) The negative tail of the 
function K c ai(<p, k c ) for five large values of k c . 



charged. Importantly, the charge distributions sharply 
peak at <p = and the angular half-width of the peak 
is about k~ . More precise, the FWHM of the peak is 
~ 2.5/ k c for k c 3> 1. Outside the close vicinities of (p = 
and 7T, the charge density U\ (tp) tends to zero for k c — > oo. 
Thus, localization of the surface charge at the sharp cor- 
ners occurs with increasing corner curvature. Remark- 
ably, the eigenfunction o\ (ip) changes its sign in between 
the charged and uncharged corners, and the correspond- 
ing zero points move towards the charged corners with 
increasing n c . This feature is not dictated by the symme- 
try properties. The higher dipolar eigenmodes 02,3.. ..(y>) 
also show localization at the sharp corners; far from them 
they possess more complicated oscillatory structure. 

Being useful to exhibit the angular localization, the 
chosen normalization is inconvenient in other respects: 
The value of the localized charge is decreasing as 1/k C7 
and the non-localized charge distributed outside the cor- 
ner areas becomes hidden. It is useful to consider the 
function K c ai((p, k c ) whose localized part is expected to 
look like the Dirac (5-function. 

Figure EDd shows in detail the negative tail of the func- 
tion K c <Ji(ip, k c ) for the same rhombus (9 a = 7r/3) and 
several large values of k c . The reduced angular inter- 
val [0, 7r/2] is fully sufficient for our analysis owing to 
the symmetry properties. All shown curves are well out- 
side the peak area, <p 3> K^ 1 . We see that after a very 
sharp initial drop, the function n c ai((p, k c ) changes its 
sign at a point (po(n c ), reaches a pronounced minimum 
at Pmin(Kc), and grows then slowly up to the second 
zero point, 7r/2. Both characteristic angles, (po(n c ) and 
Vmin^c), tend to for k c — > oo, while the minimum 
value of K c (Ti(( / 9 m in, K c ) decreases steadily with increas- 
ing k c . Far enough from the charged corner, all curves 
practically coincide with each other showing a universal 
regular behavior. Within the interval [0, 7r/2], see Fig. [3l 
the total positive (localized) charge, remains comparable 
with, but not equal to, the total negative (delocalized) 
charge. This occurs in the whole range of k c . Variation 
of the apex angle 9 p does not change this feature. 

Thus, the dipolar eigenfunctions (Tj(<p) possess compa- 
rable localized and delocalized components for k c 3> 1. 
The localized component show a complicated behavior 



in the limit k c oo. It cannot be described by the sin- 
gular function S(ip) alone. The dip in Fig. |3Jd, which 
grows steadily in the amplitude and approaches zero, 
must be responsible for an additional singular contribu- 
tion. The slowness of the transition to the limit k c — > oo 
for ai(p, k c ) correlates with that for Si(k c ). 

What is the impact of the charge localization on the 
observable characteristics of the 2D particles, such as the 
polarizability and the near-field enhancement? First, we 
determine the surface charge density <r(r) induced by an 
external light electric field of the amplitude Eo. Ex- 
panding a by the eigenfunctions o j and using the or- 
thogonality relation ((XjTji) oc Sjj/, where (...) means 
the integration along L, one can calculate the expan- 
sion coefficients [23[ . As a function of the light fre- 
quency w, the j-th coefficient is proportional to the factor 
[ej -£m(w)] _1 ; it peaks sharply at the eigenfrequency cjj 
such that e' M (ujj) = Ej provided that e^j <C \sj\. As the 
lowest value £\ (k c ) is well separated from higher eigenval- 
ues, we can restrict ourselves to a single resonant term. 
In this case we have: 
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The right-hand side does not depend on the choice of 
normalization of u \ (r) and Ti(r). With these eigenfunc- 
tions known, one calculate any linear characteristic of the 
plasmon resonance, including the dipole moment, polar- 
izability, polarization properties, extinction cross-section, 
and near fields. A similar procedure is applicable to the 
higher resonances. 

Our calculations show that the lowest mode 1 possesses 
the highest polarizability a-y for n c 3> 1. This polarizabil- 
ity exhibits no strong dependence on the corner curvature 
and remains comparable in the order of magnitude with 
the resonant polarizability of a circular cylinder of the ra- 
dius ro, |a C yi| — r o/ e M- F° r the rhombus with 6 a < n/2, 
the polarizability is anisotropic with the largest value cor- 
responding to the horizontal orientation of Eo, see Fig. [I] 
This value grows slowly with decreasing 8 a . 

The strongest impact of increasing corner curvature 
occurs for the resonant near-field enhancement. It can 
be characterized by the field enhancement factor £(<5r) = 
|E(£r)|/|E | which depends on the displacement 5r about 
the corner and also on the curvature parameter k c and 
the apex angle 6 a . We consider the optimum case when 
Eo and 5r are parallel to the long axis of a rhombus. 
Figure UK shows the product £(0)e{d versus k c exactly at 
the corner, Sr = 0, for several values of 9 a . One sees 
that £(0) is almost proportional to k c for n c ^> 1 provid- 
ing very large enhancement factors. The smaller 9 a , the 
stronger is the enhancement. The solid lines in Fig. |4Jd 
show the ratio ^ (<5t) /^(0) as a function of the horizontal 
displacement Sr for 9 a — tt/3 and three values of k c . The 
fall occurs on the scale of p c . The curves 2 and 3, plotted 
for the largest values of k c , almost coincide; they follow 
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FIG. 4. (a) The product £(0)£m versus the normalized corner 
curvature for the rhombus with different values of 6 a - (b) The 
ratio i{8r)/^{Q) versus Sr/p c for 8 a — 7r/3; the dotted line, 
(1 + 5r/p c )~ , is a coulomb fit. 



approximately the coulomb law, see the dotted line. 

Two issues are worthy of discussion - the status and 
the expected impact of this work. We suppose, we have 
filled an important gap in the knowledge of the plasmonic 
properties of nanoparticles. Two interrelated generic fea- 
tures arise from the presence of sharp corners: First, the 
positions of the resonances become predominantly depen- 
dent only on two parameters - the apex angle 6 a and the 
normalized corner curvature k c - while the shape issues 
are of minor relevance. In accordance with the Drude 
model, it causes a substantial red shift of the main dipo- 
lar resonance, X/X p w ^J2i:/9 a for n c > 10 2 , where X p 
is the plasmon wavelength. Second, strong localization 
of the plasmonic charge occurs near the sharp corners 
leading to strong near-field enhancement effects. Simul- 
taneous presence of these features allows one to excite 
the corner singularities selectively in A and control their 
strength. This can be important for such applications 
of nano-plasmonics, as nano-sensors, nano-antennas, and 
nano- heaters 0, @ ■ Nano-plasmonic nonlinear phenom- 
ena, which attract permanently increasing interest [24j . 
can greatly benefit from the presence of sharp corners; 
control of 9 a and p c is here the key issue. Last, ensem- 
bles of nanoparticles with different corners are expected 
to show very broad absorption spectra; this can be of 
interest for improvement of photovoltaic devices [25j |. 

In conclusion, we have shown that the apex angle a 
and the normalized corner curvature k c control the spec- 
trum of plasmonic resonances for 2D metal nanoparticles 
possessing sharp corners, k c 3> 1. The global form of the 
cross-section is of minor importance. The lowest resonant 
frequency of the dipolar modes is strongly separated from 
the other frequencies, and the separation grows with de- 
creasing 6 a . The corresponding plasmonic eigenfunction 
experiences progressive localization at the sharp corners 
with increasing k c . Altogether, it makes possible selec- 
tive resonant excitation of the localized states providing 
efficient near-field enhancement. The perfect-corner limit 
k c — > oc seems to be practically meaningless. 

The work is supported by the Russian Academy of Sci- 
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